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Abstract 



A promising way to reconcile naturalness with a decoupling solution to the 
SUSY flavor and CP problems is suggested by models with a radiatively driven 
inverted mass hierarchy (RIMH). The RIMH models arise naturally within the 
context of SUSY SO{10) grand unified theories. In their original form, RIMH 
models suffer from two problems: 1.) obtaining the radiative breakdown 
of electroweak symmetry, and 2.) generating the correct masses for third 
generation fermions. The first problem can be solved by the introduction of 
50(10) D-term contributions to scalar masses. We show that correct fermion 
masses can indeed be obtained, but at the cost of limiting the magnitude 
of the hierarchy that can be generated. We go on to compute predictions 
for the neutralino relic density as well as for the rate for the decay b — > 57, 
and show that these yield significant constraints on model parameter space. 
We show that only a tiny corner of model parameter space is accessible to 
Fermilab Tevatron searches, assuming an integrated luminosity of 25 fb^^. 
We also quantify the reach of the CERN LHC collider for this class of models, 
and find values of rrig ~ 1600 GeV to be accessible assuming just 10 fb^^ of 
integrated luminosity. In an Appendix, we list the two loop renormalization 
group equations for the MSSM plus right handed neutrino model that we have 
used in our analysis. 

PACS numbers: 14.80.Ly, 13.85. Qk, 11.30.Pb 
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Typeset using REVT^ 



I. INTRODUCTION 



Weak scale supersymmetry (SUSY) [|^ provides a highly motivated framework for physics 
beyond the Standard Model (SM). The hunt for the predicted supersymmetric matter is one 
of the primary tasks for collider experiments in the 21st century. The details of superparticle 
signatures are model-dependent, and may shed light on the mechanism for the communi- 
cation of supersymmetry breaking to the observable sector consisting of SM particles and 
their superpartners. 

An elegant approach is to assume supersymmetry is spontaneously broken in a "hid- 
den sector", and that supersymmetry breaking is communicated to the observable sector 
via gravitational interactions 0. Most generally, in these supergravity models, soft SUSY 
breaking (SSB) scalar mass parameters are flavor dependent [Q, and lead to violation of 
experimental limits on flavor changing (and possibly also CP violating) phenomena [Q. 
Generic predictions for CP violating, but flavor conserving, quantities such as the electric 
dipole moments of the electron or the neutron also violate experimental upper limits on 
their values. These are the well known SUSY flavor and CP problems. The proposal of uni- 
versality of scalar masses solves the SUSY flavor problem, but is not a direct consequence 
of supergravity models. Even if universality can be arranged at tree level, quantum cor- 
rections will lift the degeneracy, and lead back to problems in the flavor and CP sectors of 
the theory |P . Problems with flavor conserving, but CP violating, observables arise because 
new phases are possible in SUSY theories. These problems are ameliorated if the associated 
(combination of) phases are small, or if sparticles are sufficiently heavy. 

Over the years, considerable effort has been spent on the search for novel SUSY breaking 
schemes which are consistent with FCNC and CP violating processes. Three prominent 
classes of models include gauge- mediated SUSY breaking (GMSB) 0, anomaly-mediated 
SUSY breaking (AMSB) [|^ and gaugino-mediated SUSY breaking . Each of these models 
must introduce additional visible sector fields and/or extra dimensions, unlike the case with 
gravity-mediated SUSY breaking models. In these frameworks, the fiavor problem is solvedQ 
because particles with the same gauge quantum numbers are degenerate, but the CP problem 
may remain. 

A different solution to the SUSY flavor problem in four dimensional supergravity theories 
is decoupling 0, where the superpartner masses (for the first two generations) are raised to 
such high values that FCNC and CP violating processes are strongly suppressed. Naively, 
this solution confiicts with naturalness considerations the SUSY fiavor and CP problems 



demand multi-TeV scalar masses, while naturalness prefers sub-TeV SUSY breaking masses. 
A proposal which reconciles these requirements is to arrange for an inverted mass hierarchy 
(IMH) for scalars in the theory: multi-TeV masses for first and second generation scalars, 
but sub-TeV masses for third generation and Higgs scalars ||11|| . 

If the IMH of scalars already exists at the GUT scale (GSIMH models), then it has been 
emphasized that two-loop contributions to the renormalization group evolution of soft SUSY 
breaking scalar masses can cause a breakdown in color and/or electric charge symmetry 



^Care must be taken to ensure that the additional fields do not induce new sources of flavor 
violation. 
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unless the GUT scale top squark mass is sufficiently large |]T2|. Indeed one can generate 



viable GSIMH models only over limited portions of model parameter space. The task 



then would be to explain the origin of the peculiar choice of SSB parameters at the GUT 
scale. 

An attractive alternative has recently been suggested, and developed, in a series of papers 



1^ . The idea is to start with multi-TeV masses for all scalar fields at the GUT/Planck 
scale, and to generate the IMH radiatively (RIMH). Neglecting sub-TeV contributions to 
the renormalization group equations (RGE), and assuming Yukawa couplings and also A 
parameters remain unified at all mass scales, the one loop RGEs for the third generation 
plus Higgs scalar SSB masses can be solved explicitly. Remarkably, for simple forms of 
GUT scale boundary conditions that are consistent with SO{10) grand unification, Yukawa 
couplings rapidly drive third generation SSB and Higgs boson masses towards zero at the 
weak scale, while the first and second generation scalars (which have negligible Yukawa 
couplings) remain in the multi-TeV range. The contribution from the evolution of the gauge 
singlet right handed neutrino (RHN) superfield JSf^ that occurs in 5*0(10) unified theories 
plays a crucial role. A crunch factor S — essentially the average of first generation squared 
scalar masses compared to third generation ones — was defined. Using GUT scale Yukawa 
couplings of fcuT ~ 1 — 2, and scale choices M^r ~ 10^ — 10^^ GeV at which the RHN 
decouples from the MSSM, crunch factors of 5* ~ 50 — 800 were claimed. As an example, a 
crunch factor 5* = 400 would imply that 10 TeV ffist and second generation scalars could 
co-exist with 0.5 TeV third generation scalars at the weak scale. 

A number of problems arise if one attempts a realistic implementation of the RIMH 



model. The first, as noted also by the authors of Ref. [|T^, is that for models with large 
tan P, which is where Yukawa coupling unification occurs, it is difficult to generate a radia- 
tive breakdown in electroweak symmetry (REWSB). The second, as we will see below, is 
that for model parameters which lead to large values of S, the experimental values for third 
generation fermion masses are not obtained.^ A third problem occurs in a full implementa- 
tion of the renormalization group solution to the RIMH model. The neglected sub- TeV mass 
contributions, two-loop RGE contributions, and non-universal Yukawa couplings and A pa- 



rameters below the GUT scale all serve to perturb the simple solutions found in Ref. |0 



generally resulting in a considerable diminution of the crunch factors that can be obtained. 



It has been shown [|1^] that the problem of implementing REWSB in Yukawa unified 
5*0(10) models can be solved by including allowed D-term contributions to scalar field 
masses: these occur when spontaneous gauge symmetry breaking leads to a reduction in 
rank of the gauge group [1^. The form of the 5*0(10) Z^-terms allows for a split in the GUT 
scale Higgs masses so that frifj^ < which, in turn, facilitates REWSB. The D-terms also 
affect other scalar masses, resulting in a sparticle mass spectrum where the bottom squark 
bi is frequently the lightest of all the matter scalars, and left- sleptons are lighter than right- 
sleptons, as opposed to the situation in the mSUGRA model, with GUT scale universality 

In a recent paper |T^, it was shown that 50(10) D-terms could rescue REWSB into the 



^Yukawa couplings, and hence fermion masses, for the first two generations, and CKM matrix 



parameters, are neglected in Ref IQ, in Ref. IT^], and in this paper. 
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RIMH framework as well, so that sparticle masses could be explicitly calculated. In Ref. ||17||, 



a bottom up approach was used to calculate the values of GUT scale Yukawa couplings, and 



the GUT scale SSB boundary conditions identified in Ref. ||T^ were adopted. The sub-TeV 
contributions to RGEs were also included. The crunch factors obtained, however, were in 
the range S* ~ 2 — 7, much smaller than the values found previously |T^ . 



In this paper, we perform a more detailed analysis of RIMH models. We have two main 
goals. First, we seek a realistic implementation of RIMH models, keeping track of sub-TeV 
and two-loop contributions to RG evolution, requiring appropriate values of third generation 
fermion masses, and requiring a radiative breakdown in electroweak gauge symmetry. We 
find it possible to construct realistic RIMH models, but with crunch factors 5* ^ 10. Thus, 
we conclude that realistic RIMH models give only a partial solution to the SUSY flavor and 
CP problems. The RIMH model, coupled with a limited degree of non-universality of soft 
breaking scalar masses or alignment of fermion and sfermion matrices may be enough to 
obtain a complete solution. Once realistic SUSY particle mass spectra can be calculated 
within the RIMH model context, our second goal is to examine various phenomenological 
implications of these models, including the neutralino relic density, b s'j decay rates, and 
expectations for collider searches. 

Toward these ends, we upgrade the supersymmetric RG solution used in ISAJET to 
include the complete set of two-loop RGEs for both the MSSM and MSSM-I-RHN models. 
As first shown in Ref. [Q, multi-TeV two loop contributions can potentially be as large as 



sub-TeV leading order contributions to the RGEs. In Sec. Ha, we first re-analyze RIMH 
models using a bottom-up approach to the SUSY RG solution. The bottom-up approach 
guarantees solutions with correct third generation fermion masses. Allowing more flexibility 



in the value of Mj^ allows somewhat greater values (compared to Ref. [|T^]) of S* ~ 9 to be 
obtained. These S values are far below the values of S claimed in Ref. |T^. To understand 
this, in Sec. 116, we repeat much of our analysis using a top-down RGE approach. In this 
case, large GUT scale Yukawa couplings can be used. We find then that a significantly 
larger IMH mass gap, with S* ~ 35 can be obtained. However, in this case the third 
generation fermion masses do not match their measured values. In particular, the fixed 
point prediction of the value of rrit is quite far from its experimental value. The remainder 
of the paper focuses upon the phenomenological implications of the model. In Sec. Ill, we 
attempt to systematize the parameter space explored in Sec. II, and show plots of sparticle 
masses versus model parameters. In Sec. IV, we present our predictions for the cosmological 
neutralino relic density, while those for the branching ratio of the 6 — 57 decay are shown 
in Sec. V. Finally, in Sec. VI, we show the reach of the Fermilab Tevatron and CERN 
LHC for RIMH models. In spite of the decoupling solution with multi-TeV first and second 
generation scalars, a considerable reach is obtained at the CERN LHC, corresponding to 
values of rrig ~ 1600 — 1700 GeV, assuming just 10 fb~^ of integrated luminosity. We end 
with a summary of our results along with some conclusions in Sec. VII. In an appendix, we 
list the two-loop RGEs for the MSSM and MSSM-I-RHN models used in our analysis. 

II. RIMH MODEL 



We begin by briefly reviewing the analysis of Bagger et al. \14 . Within the framework 



of the MSSM-I-RHN model (which is assumed to be valid up to the GUT scale), the 1-loop 
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RGEs for the third generation and Higgs scalar masses and trilinear terms can be written 

as 



jrn^ = -rNm2, (2.1) 

where = {m'jj^,mfj,mQ,mj^,m'jj^,m'j^,m'^^,m'lj, A^)"^ , t is the natural log of the energy 
scale and Y = -j^. In Eq. (pTll), sub-TeV scale terms such as gaugino masses are neglected. 



Inspired by 50(10) unification, / here is the unified Yukawa coupling, and A is the unified 
trilinear SSB coupling and N is a matrix of numerical coefficients obtained from the one- 
loop RGEs. The qualitative behaviour of the weak scale SSB masses can be obtained by 
neglecting the splitting between the third generation Yukawa couplings and A parameters 
that arises from renormalization group running. The evolution of the SSB masses is then 
given by. 



r /■*/ 

m^(t = tf) = Cjejexp —Xi J Y dt 



(2.2) 



where Aj is the ith eigenvalue corresponding to eigenvector ej of the matrix N, and tf is the 
natural log of the weak scale. 

The eigenvalues of the matrix N are non-negative, with the largest eigenvalue Ai = 
16. Components of the initial conditions along this direction are strongly suppressed by 
renormalization group evolution to the weak scale: the corresponding eigenvector suggests 
that the SSB masses should be related as 

Ami = ^"^u = 4:ml = Ami = 4m| = Amj, = 1m\^ = 1m\^ = Al, (2.3) 

at the GUT scale. This condition is consistent with minimal 5*0(10) unification. Plugging 
the eigenvector ei and eigenvalue Ai into (|2.2|) , it is found that third generation SSB masses 
are exponentially "crunched" to zero at the weak scale tf. First and second generation 
scalars have only tiny Yukawa couplings, and hence, their masses are essentially unchanged 
by the evolution between the GUT and weak scales. This is the basic mechanism for the 
radiative generation of the IMH. The singlet neutrino superfield N!^, which is automatically 
present in 50(10) SUSY GUT models is critical for this mechanism to work. 

In our examination of the RIMH model, we adopt several assumptions motivated by the 
work of Bagger et al, Ref. |T^. First, we assume the existence of an 50(10) SUSY GUT 



model above the scale Q = Mqut, where Mqut is defined as the scale at which the gauge 
couplings gi and g2 unify. The 50(10) gauge symmetry breaks to 5[/(3)c x SU{2)l x f/(l)y 
at Mguti so that the MSSM-I-RHN model is the effective theory immediately below the GUT 
scale. The matter superfields of the model include the usual superfields of the MSSM plus 
a gauge singlet neutrino superfield: 




Q^ = \l \ . u = ri]. u::, ei ni 



c 



where z = 1, 2, 3 corresponds to the various generations. In addition, the Higgs multiplets 
are given by 
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i/„(2)^(|).a„diJ.(2)^(|). 

The superpotential is given by 

/ = fussM + fue^.VHiN' + ^M^N^N', (2.4) 

where 

Imssm = fJ-H^Hdi + fu^ijQ'^HlJJ'^ + fdQ^HdiD^ + feL^HcuE'^. 

eij is the completely antisymmetric SU{2) tensor with ei2 = 1, and we neglect terms in- 
volving Yukawa couplings of the first and second generation. Here, and in the following, we 
show terms involving only third generation right handed neutrino and sneutrino fields, and 
suppress the generation index on these. 

The soft SUSY breaking terms can be written as 



^ — C-MSSM — ^\p\^r\^ + 



AJ^eijVHiuli + ^B^MNi>l + h.c. 



(2.5) 



where 



jCmssm = - E ^l\^r\^ M^KX^ + [BiiHlHdi + h. 

r ^ \ 



c. 



The parameters A^, By and rriy^ are assumed to be comparable to the weak scale. Below 
the scale Q = M^, but above the weak scale, the RHN superfield is integrated out, and the 
MSSM is the effective theory. In our renormalization group analysis, we use the full set of 
two loop renormalization group equations for the MSSM and the MSSM+RHN models, as 
given in the Appendix. These RGEs have been encoded into the event generator ISAJET 
7.51 [jl8i, which we use for this analysis. 



A. Bottom-up approach 

To generate supersymmetric spectra within the framework just described, we adopt the 
bottom-up approach inherent in ISAJET. This ensures that the matter fermion masses are 
as given by experimental measurements. Beginning with weak scale central values of the 
three SM gauge couplings and three third generation Yukawa couplings, we evolve upwards 
in energy until gi = g2, where the GUT scale is identified. In keeping with the spirit of 
5*0(10) unification, we then take the various GUT scale SSB scalar masses to be 

ml = ml = ml = mle + 
ml = m\ = m\^ — 3M|, 
= m^g + 5M^ 
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where mie is a common mass of all scalars in the 16- dimensional spinorial representation of 
5*0(10), and mio is the mass of the Higgs scalars, which are assumed to reside in a single 
ten-dimensional fundamental representation of SO{10). The parameter M^i reflects the 
magnitude of the D-teim contribution to scalar masses induced by the breakdown of 5*0(10); 
its magnitude is taken to be a free parameter. We also assume the gaugino masses and 
trilinear SSB terms unify to mi/2 and Aq at Mqut, respectively. We do not impose Yukawa 
coupling unification at Mqut, although as we shall see, the largest IMHs are generated 
when Yukawa couplings most nearly unify. Beginning at Mguti we evolve the SSB masses, 
gauge and Yukawa couplings down in energy to a scale Q = Mjy using the full two loop 
MSSM+RHN RGEs given in the Appendix of this paper. From Q = Mj^ to M^^ak, we 
evolve according to two loop MSSM RGEs. At M^eak, radiative electroweak symmetry 
breaking is imposed, using the renormalization group improved one-loop scalar potential, 
and SUSY loop corrections are included for the relation between matter fermion masses 
and the corresponding Yukawa couplings [0. An iterative procedure is used in ISAJET to 
obtain the superparticle mass spectrum. 

The model is completely specified by the parameter set, 

mi6. Ml,, mi/2. Mat, tan/3, sign{fi), 

where mio and Ao are determined in terms of mig by the boundary condition ( |2.3| ). We 
adopt pole masses = 1.777 GeV, nib = 4.9 GeV and rrit = 175 GeV. We generate random 
samples of model parameters with the following ranges: 

500 < mi6 < 10000 GeV, 

< Ml < mlj3, 
100 < mi/2 < 3000 GeV 
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10^ <Mn < 10^^ GeV, 
10 < tan/3 < 55, 
/i > or /i < 0. 

Our first results are shown in Fig. ||. Here, we show the "crunch" factor 5" is defined as, 
3 (m? + m J + m? + m j ) + m? + m? + m? 

g _ ^ d-L dfi' <^L ' Efi ' '^e (2 6) 

3(m| + m?^ + m| + m?^) + m\ + + m?^ ' 
for solutions to the superparticle mass spectrum consistent with REWSB. We take /i < 



smce 



and < in this figure. Our definition of 5* differs slightly from the one in Ref. [0 
we are able to use mass eigenvalues in the calculation. In frame a), we show 5* versus mig. 
We first see that the largest values of 5 obtained are close to ~ 7. In dedicated searches, 
values of 5 as high as 9 have been generated. The highest 5 values occur for mig ~ 2000 
GeV. For higher values of mig, two loop RGE can become important, thus perturbing even 



more the idealized solution of Ref. [0], and possibly leading to tachyonic mass parameters 
|l^,13]. In frame &), we show 5* versus tan/?. Here, it is easy to see that the greatest IMH 
develops for very large values of tan/3 ~ 50 where Yukawa unification can occur. Frame 
c) shows 5* versus a ratio R indicative of the degree of Yukawa coupling unification. We 
define Rth = rnax{ft/ fb, ft/ ft), with similar definitions for Rtr and Rrb, and where ft, fb 
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and fr are the third generation Yukawa couphngs evaluated at Q = Mqut- Then, R is the 
maximum of i?,-?,, Rtb and Rtr- The largest S values occur at values of i? ~ 1, showing 
that the largest hierarchy is obtained for models with the best Yukawa coupling unification. 
Finally, in frame (i), we plot model solutions versus the superpotential neutrino mass M^v- 
In this case, there is a mild correlation to obtain higher S values at small to intermediate 
values of Mat, although a significant IMH can be generated at any of the values shown. 

In Fig. 0, we plot in a) S versus the ratio mi/2/mi6. In this case, we see 5* reaching 
a maximum for mi/2/''^i6 ~ 0.2. If mi/2 is comparable to mig, then the effect of gaugino 
masses will be important, and will destroy the radiative IMH solution. If mi/2 is too low, 
then would-be large 5* models are excluded because we require that the lightest SUSY 
particle (LSP) should be a neutralino: in these scenarios, for a fixed value of mig, third 
generation sfermion masses decrease faster with mi/2 than does Tn^^. In frame fe), we show 

S versus the ratio Mu/rriiQ. The distribution is bounded from above by l/\/3; higher values 
imply tachyonic sbottom masses at the GUT scale. We see that 5* reaches a maximum at 
Md/ttliq ~ 0.2. Models with a high value of tan/? required by Yukawa unification need a 
non-zero value of to achieve a consistent REWSB mechanism; models with Md ~ are 
allowed in our plot, but will typically have smaller tan/3 values, and hence a smaller degree 
of Yukawa unification. 

In Figs. ^ and ^, we show similar plots, but for > 0. Many of the results are 
quahtatively the same as in Figs. [| and ^. An exception occurs in the S versus R plot, 
where we see that values of i? ~ 1 are not allowed. This is a reflection of the well-known 
lack of Yukawa coupling unification for /i > in a bottom-up approach. We have also made 
similar plots for Aq > 0. In general, these models lead to smaller values of S, and so are less 
promising than the Aq < models. We do not show these for reasons of brevity. 

Our conclusion is that the maximum crunch factors obtainable in realistic RIMH models 
are typically S* ~ 7 — 9. While these crunch factors are considerably smaller than values 
found in Ref. [l^, they can still lead to models with multi-TeV first and second generation 
scalar masses, which is sufficient to solve several, but not all, of the problems associated 
with FCNCs and CP violating processes. 

In Table |, we show two sample RIMH model solutions, taking mig = 2500 GeV, mi/2 = 
650 GeV, tan/? = 50, and = 500 GeV. The first solution, for /i < 0, has Yukawa 
coupling unification to ±1.5%, with GUT scale Yukawa couplings at / ~ 0.54. The bi 
squark is the lightest of the third generation scalars, with = 446.3 GeV, while first 
and second generation scalars have masses around 2400 — 2800 GeV. Note the relatively 
heavy light Higgs scalar rrih = 129.1 GeV, which may make detection at Fermilab Tevatron 
experiments difficult. The tau sleptons have mass just beyond 1 TeV, while the lighter top 
squark has m^-^ = 858.8 GeV. 

The second case study has all the same parameters, but the sign of /i is flipped: /i > 0. 
In this case, SUSY loop corrections to Yukawa couplings make Yukawa unification difficult, 
with ftiMcur) = -495, while fbiMcuT) = -305. The bottom squarks are all beyond 1 TeV, 
but = 742.0 GeV, and the light tau slepton is the lightest third generation scalar with 
= 661.2 GeV. For both cases, first and second generation scalars are heavy enough to 
suppress many flavor changing and CP- violating effects, but not all of them. Such a model 
would need to be coupled to some limited universality to solve, for instance, the SUSY flavor 
and CP problems in the kaon sector. 
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B. Top-down approach 



In the bottom-up approach of the previous subsection, starting with central values for 
the weak scale Yukawa couplings, we only were able to obtain GUT scale Yukawa couplings 
of order / ~ 0.5. These values of GUT scale Yukawa couplings are significantly lower than 
the choices used in Bagger et al, where /{Mgut) ~ 1.0 — 2.4. Extrapolating their results 
to lower values of GUT scale Yukawa couplings indicates that S ^ 25 — 30, only about 
a factor 3 larger than what we find. We attribute this difference to the additional effects 
such as inclusion of sub-TeV gaugino mass terms and two loop effects in the RGEs, as well 
as deviation from the ideal boundary condition due to the D-teim contributions to scalar 
masses, all of which would result in a smaller crunch factor. 

To examine why we were confined to small values of the unified Yukawa coupling, we 
modified our computer code to adopt a top-down RGE solution in ISAJET. In this case, we 
could use arbitrary values of a unified GUT scale Yukawa coupling as inputs, and compute 
the masses of third generation fermion (as well as their superpartners). Our results for fi < 
are exhibited in Fig. |^. In frame a), we show the pole mass rrit versus the value of f{MGUT)- 



The well-known Yukawa coupling fixed point behavior is evident in this plot, where a 
large range of values of /(Mcc/t) = 1 — 2 leads to rrit ~ 200 — 210 GeV. We see, however, 
that values of fi^McuT) — 0.5 are necessary in order to obtain rrit = 175 GeV, and further, 
that for /{Mgut) > 1, ""^t > 190 GeV, which is conclusively excluded by experimental 
measurements.^ 

In frame b), we show the crunch factor 5* as a function of f{MGUT)- Taking fiMcuT) — 
0.5, which leads to the measured central value of the top quark mass, yields 5^5. Taking 
much larger values of fiMcuT) can lead to a much greater IMH, with ~ 25 — 35 being 
possible for f^Mcur) — 1-2, to be compared with 5* ^ 100 in Ref. [|I4|- However, such 
models will in general lead to a value of the top quark mass significantly above its measured 
value. For yet larger values of /{Mgut), electroweak symmetry is not correctly broken, and 
these large 5* solutions are no longer obtained: for solutions with tan /3 near the maximum 
of its allowed range, we get /z^ < when /{Mgut) is increased; for other cases, we find that 
third generation mass parameters become tachyonic.0 Frame c) shows S versus the ratio 



^The reader might worry that we have used 1-loop relations for the connection between pole and 
running top masses, or alternatively, that these might be scheme dependent. We have checked that 



in the MS scheme, the difference between the 1- and 3- loop [21| running top masses is < 2.2 GeV, 
and further, that the difference between the running masses in the AIS and DR schemes [^] is 
< 2 GeV. As stated previously, we have included 1-loop SUSY radiative corrections |l^ to the 
relation between the Yukawa couplings and the corresponding fermion masses, and found these to 
be significant. 

^The precise value of /{Mgut) for which this occurs is somewhat sensitive to the procedure that we 
adopt for numerically integrating the RGEs. We take all sfermion masses to be mig during the first 
iteration of the RGEs. Especially in the case of the IMH model, this grossly overestimates the scale 
Q at which the Higgs potential parameters are frozen in ISAJET. In some cases though this results 
in tachyonic third generation masses at this stage, and the parameter point is rejected. While this 
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Md/ttliq in the top-down approach. The largest 5* values are obtained for smaller values of 



the SO{10) D-teim. Larger D-terms not only perturb the boundary conditions of Ref. [14 



but can lead to tachyonic bottom squark masses. Finally, in frame d), we show S versus 
the ratio ^«i/2/'"^i6- In this case, there is a tight correlation, with the largest S values being 
obtained for small values of mi/2- Larger mi/2 values again will increasingly perturb the 
simple evolution described in Ref. [1^, which was obtained assuming that gaugino mass 
could be neglected. We have also analyzed the corresponding situation for /i > 0. The 
results, which are shown in Fig. ^, frames a)-d) are very similar to those in Fig. ^ just 
discussed. 

We illustrate the spectrum of a model with a large unified GUT scale Yukawa coupling 
which yields S = 33.8 in Table ||. Here, we take f^Mcur) = 1-28, with other parameters 
as shown. In this case, first generation masses of ~ 12.5 TeV are generated, while third 
generation scalar masses are typically 1 — 3 TeV. The spectrum is meant for illustrative 
purposes only: not only is the top quark mass for this case mj = 204.8 GeV, but third 
generation scalars are too heavy to claim that the weak scale is really stabilized. 



III. SPARTICLE MASSES IN THE RIMH MODEL 

In this section, we attempt to systematize the parameter space and map out sparticle 
masses in preparation for a phenomenological analysis of the RIMH model. /^From Figs. |l| 
and D, we see that the largest crunch factors are obtained for large tan/3 values of ~ 50, 
where Yukawa coupling unification occurs with the best precision. In this region of parameter 
space, D-teim contributions to scalar masses are needed to gain consistency with REWSB, 
and in fact we see from Figs. and|^ that values of ~ 0.2mi6 lead to the largest crunch 
factors. 

Motivated by these considerations, we show in Fig. the allowed regions of the 
rriiQ vs. mi/2 parameter plane, taking tan/? = 50, Mo = 0.2mi6, = 1 x 10'^ GeV 
and /i < 0.0 Throughout the rest of the paper we take Aq < 0, since this generally allows 
a bigger IMH. For the lattice of points in this plane we have used ISAJET to obtain the 
sparticle masses, subject of course to several theoretical constraints. Points with triangles 
are excluded because the LSP is not a neutralino, points with open circles are excluded be- 
cause they do not lead to REWSB (/^^ < 0), while those denoted by dots lead to tachyonic 



is indeed an artifact of our calculational procedure, we do expect tachyonic third genration mass 
parameters when the Yukawa coupUng becomes sufficiently large. Thus, while the precise value of 
f{MGUT) beyond which there are no large S solutions may be somewhat different from that shown 
in the figure, we believe that the figure is qualitatively correct. Since the value of f{MQUT) at 
which the large S solutions really cuts off is only a matter of academic interest, and has nothing 
to do with the real phenomenology of the model, we have not pursued this any further. 

^Our choice of means that the simple see-saw mechanism (with a generation-independent 
RHN mass) will not give neutrino masses in agreement with the Super Kamiokande data. There 
must be some additional physics in the neutrino sector which, we presume, will not significantly 
affect sparticle phenomenology. 
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scalar mass parameters. We show contours of regions where S > 3 and 5* > 5, indicating 
that S increases as one proceeds towards the lower-right region of the plot. Values of 5* as 
high as 8 can be obtained adjacent to the lower right excluded region. Contours of rrig = 2 
TeV and m^^ = 3 TeV are also shown, with particle labels on the contour side with lower 
mass values. We also show contours of m^-^ and mj^ of 1 TeV, and fi = —1 TeV. The regions 
below the ti and bi contours are most natural, while the region above the ul contour best 
suppresses FCNC and CP violating loop effects. Finally, a contour shows the region where 
R < 1.05, corresponding to Yukawa coupling unification within ±2.5%. 

To obtain an idea of how some of these quantities change with M^, we plot a) S, b) R, 
and d) m^^ versus mig for mi/2 = 0.25mi6 in Fig. ||for several values of indicated 



c m 



on the figure. The curves cut off when they run into the excluded regions discussed in Fig. |^. 
/,From a), we see that, for a given value of mie, the largest S values are obtained for the 
smallest values of Mn- From b), the smallest Mn values also lead to models with the best 
Yukawa coupling unification. In c), it is seen that low values of miQ and can lead to very 
small values of , some of which might even be accessible to direct bibi search experiments 
at the Fermilab Tevatron [^. In d), the lightest top squark masses are generated for small 
mi6 but larger values of M^, so that some parts of model parameter space may be accessible 
to direct titi search experiments at Fermilab Run 2 p4|. 

In Fig. ^ we show various sparticle and Higgs boson masses and the fi parameter versus 
rriiQ for the same parameters as in Fig. |^, but with Mjv = 1 x 10^ GeV.Q In frame a), a 
significant mass gap is seen between first and third generation scalars. Although the value 
of S decreases as mig increases, the mass gap between the two generations is increasing. 
For the first generation sleptons, we see that m^^ < mgj^, as opposed to the situation in 
mSUGRA models with universality of scalar masses. This can be a result of both D-term 



contributions to scalar masses, and the effect of two-loop contributions to the RGBs |13|. In 



frame 6), we see that while rrig and are increasing with mig, the value of is decreasing, 
and in fact where /i — > forms the upper limit of the mie parameter space. This impacts 
upon the nature of the two lighter neutralinos and the light chargino, which are gaugino- 
like for small mie, but become increasingly higgsino-like for larger values of mig, and has 



implications for the cascade decay patterns of gluinos, which should be accessible ||26 
at the LHC. 



In Fig. |T0|, we show again the mig vs. mi/2 plane, as in Fig. 0, but for /i > 0. In this 
case, much more of the parameter plane is excluded. The measure of Yukawa unification R 
varies from 1.6—1.9 across the allowed region, and a slightly reduced value of crunch factor 



S is obtained compared to the negative fi case. In Fig. |TT|, we show the corresponding 
values of S, R, m^^ and m^-^ versus mig, with the same parameters as in Fig. ||, except 
mi/2 = 0.22mi6, and /i > 0. In this case, the largest S values are obtained for the largest 
values of M^, which is also where better Yukawa unification is achieved. The bi mass here is 
always rather high, and not accessible to Tevatron searches. However, the ti is the lightest 



^ The calculation of /i near the edge of the excluded region (where \fi\ 0) is numerically delicate, 
so we have increased the number of Runge-Kutta steps to 1000 and convergence limits to 0.002, 
beyond their default ISAJET values, for these plots. 
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of the scalars, and may be accessible to Tevatron searches at Run 2 if mig is small enough. 
Its mass is almost independent of the value of Mjv- Finally, in Fig. 0, we show again 
various sparticle masses for the same parameters as in Fig. |Tl|, but with mi/2 = 0.22mi6 
GeV. Once again, we see that the light chargino and the two lighter neutralinos go from 
being gaugino-like to higgsino-like as mig increases. 

Before moving on to the phenomenology, we remark that several authors have suggested 
that since n enters the tree level expression that determines Mz via the REWSB constraint, 
H ^ 300 GeV might provide a more valid determination of naturalness than the third 
generation scalar masses. Adopting this criterion, we would then conclude that in the 
region of the parameter plane in Fig. ^ or Fig. [1^ close to the open circles (where yU^ 0) 
the theory would be technically natural. We note, however, that n generically tends to be 
large, and very sharply dives to zero over a limited range of parameters. This means that 
a small change in, for instance, mig could lead to quite a different value of fi in this region. 
For this reason we do not adopt as a naturalness measure in our analysis. 



IV. NEUTRALINO RELIC DENSITY 

The calculation of the cosmological relic density of neutralinos ^^_^h'^ provides an im- 
portant constraint on supersymmetric models with i?-parity conservation. ^^Ji^ > 0.02 to 
account at least for galactic rotation curves. Analysis of the cosmic microwave background 
suggests a total energy density of the universe = 1.0 ± 0.2. With the Hubble parameter 
estimated to be about 65 km/sec/Mpc, we have Qh'^ ~ 0.42. In a universe with 5% each of 
baryonic matter and massive neutrinos, and a cosmological constant with Q/i^h"^ ~ 0.25, we 
would expect ^^_^h'^ ~ 0.13. As an extreme, assuming no contribution to dark matter from 
the cosmological constant, we would obtain ^j;Ji^ ~ 0.38. 

Our estimate of the relic density of neutralinos in the RIMH model follows the calcula- 
tional procedure outlined in Ref. [^. Briefly, we evaluate all tree level neutralino annihila- 
tion diagrams exactly as helicity amplitudes. We then calculate the neutralino annihilation 
cross section, and compute the thermally averaged cross section times velocity using the 



fully relativistic formulae of Gondolo and Gelmini ^8[. Once the freeze-out temperature 
is obtained via an iterative solution, we can straightforwardly obtain the neutralino relic 
density VL^Ji^. Our program takes special care to integrate properly over any Breit-Wigner 
poles in s-channel annihilation diagrams [^, an important feature to reliably estimate the 



relic density at large values of the parameter tan (3. We neglect co-annihilation effects, which 
can be important when the Zi is nearly mass degenerate with any of the other sparticles 



Our results for the neutralino relic density are shown in Fig. |T3|. We plot ^'^_^h versus 



for both signs of /i. Here, we take tan/5 = 50 and Md = 0.2miQ. For /i < 0, we 
take mi/2 = 0.25mi6, while for /i > 0, we take mi/2 = 0.22mi6, i.e. we calculate along a 
diagonal strip in each of Figs. |^ and 0, in the region where S is large. For /i < 0, values of 
mi6 ^ 2200 GeV are in the excluded region as can be seen from Fig. |^. 

Usually, in the mSUGRA model with large values for scalars masses at the GUT scale, the 
relic density is too large as a result of scalar mass suppression of the neutralino annihilation 



cross section. In fact we see from Fig. 13 that significant ranges of mig can give rise to 



either a too small or a too large neutralino relic density. In the case that /i < 0, only a very 
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narrow region with mig ~ 4400 — 4600 GeV gives 0.13 ^ ^zj^"^ ~ 0.38, the favoured range. 
In this region, Zi is a mixture of higgsino and the bino, and annihilates efficiently via Z 
exchange, until for very large mie the LSP is dominantly a higgsino, and the relic density 
becomes too small.[] Conversely, for mig smaller than 4400 GeV, the LSP is dominantly 
bino-like (see Fig. |^&), and the relic density rapidly increases because scalars are heavy and 
the LSP couplings to gauge particles are suppressed. For yet smaller values of mie, the 
curve turns over and the relic density again drops, mainly because of efficient annihilation 
into bottom pairs via the exchange of a light hi (see Fig. |^a) in the t- and u-channels. In 
the /i > case, our relic density computation favours the region around mig = 3600 GeV. 
For larger values of mie, the relic density again becomes too low for the same reason as for 
the /i < case. For smaller values of mig, the relic density rapidly becomes too large, but 
again turns over for mig < 3000 GeV when the bottom and top squarks are light enough 
to allow significant LSP annihilation to hb and ti pairs, respectively. Indeed, the additional 
turnover in the curve just below mie = 2000 GeV is exactly where Zi annihilation into 
top pairs becomes possible. In summary, the relic density constraints favour a relatively 
narrow range of mig where the LSP is in the bino-higgsino transition region. Furthermore, 
significant portions of parameter space in the RIMH model, yield VLh"^ > 1, and so lead 
to too young a universe, and are excluded. However, viable parameter space regions can 
remain when certain squark masses are sufficiently light, or where the Zi is in the transition 
region between being bino-like or higgsino-like. 

We remark that in the parameter space region with a higgsino-like LSP, it is important 
to note that rates for neutralino-nucleon scattering are also enhanced and this en- 
hancement is compounded at large tan /3. Thus, in the higgsino-like LSP region (where the 
cosmogical dark matter must have a different origin), direct dark matter search experiments 
can hope to find the first evidence of supersymmetry by detection of neutralino-nucleon 
scattering. 



V. THE DECAY & ^ 57 IN THE RIMH FRAMEWORK 

It is well known that the decay 6 57 provides a stringent test of physics beyond the 
SM. The reason is that even within the SM framework, this decay can only occur at the 
1-loop level, via a loop involving a top quark and a W boson, with the photon attached 
to either one, and so, the SM amplitude for this is often comparable in magnitude to the 
corresponding amplitude from new physics. The CLEO experiment still provides the 
best determination of its branching fraction, 

BR{b 57) = (3.15 ± 0.35 ± 0.32 ± 0.26) x 10"^ (5.1) 

and restricts the branching ratio at 95% CL to be 

2 X 10"^ < BR{b 57) < 4.5 x 10-^ (5.2) 



' Neutralino co-annihilation will further diminish the relic density in this region. This may slightly 
shift the allowed range of mig from that in the figure. 
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while the precision, 



BR{b 57) = (3.34 ± 0.5 ± 0.36 ± 0.27) x 10-^ (5.3) 



from the BELLE experiment |^T| is not much less. These results are in excellent agreement 
with the SM prediction, which at NLL accuracy |3^ is. 



Si?(6 ^ S7) = (3.28 ± 0.33) X 10-^ (5.4) 

The calculation of the width for 6 — S7 decay proceeds by calculating the loop interaction 
for 6 — i> 57 within any particular model framework, e.g. the MSSM, at some high mass scale 
Q ~ Mw, and then matching to an effective theory Hamiltonian given by 

Heff = ~^V,,VZJ2C^{QmQ), (5.5) 
v2 i=i 

where the Ci{Q) are Wilson coefficients evaluated at scale Q, and the Oj are a complete set 
of operators relevant for the process b 37. We have included all order QCD corrections 
via renormalization group resummation of the leading logs which arise from the disparity 
of the new physics scale and the scale rrif, relevant to the decay. Our procedure has been 



described in detail elsewhere |^3|], and we will not repeat it here. 

However, to enable the reader to get some feel for our results, we mention that within 
the MSSM, the Wilson coefficients C7 and Cg receive additional (1-loop) contributions tH~ , 
qiWj, Qig and (jiZj loops, but the first two, by far, dominate at least within the mSUGRA 
framework. We should think of these coefficients as being renormalized at a scale close to 
the scale of sparticle masses. As mentioned above, the coefficients have then to be evolved 
down to Q = rrih. The amplitude for the decay is proportional to Cj{mi^ which, in turn, 
depends on Ci{Mw) {i = 2, 7, 8). 

We now turn to the results of our computation of BR{b — *• S7) within the RIMH frame- 
work. We have computed this branching fraction for the same two "large S" slices as in 



Fig. |T3|. Our result is shown in Fig. |14|. We see that for the slice with negative /i, the branch- 



ing ratio that we obtain is considerably higher than the 95%CL upper limit ([5. 2]), except 



for the largest values of tuiq in this figure. If we recall that large S models with negative fx 
tend to give better Yukawa coupling unification, the reader will immediately recognize that 
the confiict between the experimental value and our prediction for the 6 — 57 decay has 
been noted before: models with the sign of fi that allows Yukawa coupling unification tend 
to give large values for BR{b 57) 0,0. We note though that for ttiiq ~ 4600 GeV, 
the value favored by our analysis of the relic density, the branching ratio is just about 25% 
larger than the CLEO upper limit. It is possible to imagine additional effects {e.g. squark 
mixings or SUSY phases) that could serve to reduce it to an acceptable level. 

For the slice with > 0, we see that the RIMH model is in excellent agreement with 
experiment for mig ~ 1800 — 2000 GeV. Unfortunately, this region of mig seems to be 



disallowed by our computation of the relic density from neutralinos in Fig. |I3|. For smaller 
values of mig, the branching ratio rises rapidly to even beyond its maximum for the negative 
fi case, while for larger values of /i, it decreases to well below the experimental lower bound 
from (|5.2|). An especially intriguing feature is the sharp dip near mig = 3700 GeV. 
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To gain some understanding of the behaviour seen in Fig. |TJ, we have examined the 
dominant contributions to the Wilson coefficients C-j and Cg (as mentioned C2 does not get 
any SUSY correction so that C2{M\y) ~ 1) as a function of mig for the same parameter 



space shces. We show our results for the negative and positive /x cases in Figs. |T5| and 
Figs. |1^, respectively. In frame a), we show our computation of C7 while the corresponding 
computation for Cg is shown in frame h). The dashed lines show the dominant individ- 
ual components, while the solid lines labelled Cj{Mw) and Cs{Mw) show the sum of all 
individual contributions. Finally, we have also shown the result after evolution of these 
coefficients down to the weak scale. As mentioned, Cj{mi,) is the important quantity for the 
computation of the rate for the decay b sj, and it is for this computation that C2^s{Q) 
are needed. 

The following points are worth noting. 

1. Contributions from gluino and neutralino loops are negligible and are not shown. 

2. The contribution from the H"^ loop is small, and only weakly dependent on mie, 
because mfj± (which is close to 171^) is always large. The SM contribution is, of 
course, independent of mig. Essentially all the dependence on mig comes from the 
tiWj contributions, which tend to grow as mig gets close to its lower end, primarily 
because the top squarks tend to be relatively light, as can be seen from Fig. ^ and 



Fig. 12 



For /i < 0, the SUSY contributions interfere constructively with the SM contribution. 
Since the latter, by itself, gives good agreement with the experimental decay rate, it 
is not surprising that the SUSY model prediction for the rate is too large for this sign 
of /i. 

The SUSY amplitudes, which for large tan/3 are proportional to Atfitanf], reverse 
their sign if /i > 0, and interfere destructively with the SM amplitude. Indeed it is this 
that makes it possible for the RIMH model to be in agreement with the experimental 
value for BR{b 37) . For this to be possible, the RIMH model amplitude (aside from 
the QCD evolution effect) has to have about the same magnitude (but opposite sign) 
as the amplitude in the SM. 

We see that the chargino contribution, and thus the total amplitude, varies much more 
for the positive /i case than for the case with negative /x. Presumably, this is because 
top squarks can get much lighter when /i > as can be seen from Fig. This also 



explains why the range for the prediction in Fig. |I^ is much wider for the positive 
case. The sharp dip in the prediction of the branching ratio occurs when Cj{mb) goes 
through zero. 

The evolution of C-j between Q = Mw and Q = is much more sizeable for the 
fi > case than for /i < 0. Indeed we see that this effect reduces T{b — > 57), which 
should be proportional to IC7P, quite substantially when /i > 0. 
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VI. REACH OF FERMILAB TEVATRON AND CERN LHC FOR RIMH MODELS 



Supersymmetric models with an IMH spectrum for the superparticles present unique 
problems and opportunities for detection at collider experiments. Since scalars of the first 
two generations have multi-TeV masses, they effectively decouple, and will be produced 
(if at all) with only tiny cross sections. Charginos, neutralinos and gluinos may also be 
too heavy, so that search experiments may be forced to focus on direct detection of third 
generation scalars. In addition, third generation sleptons are generally difficult to detect 
at hadron colliders, so that searches may have to focus on direct production of top and/or 
bottom squarks. 

If /i < 0, then from Fig. |^ we see that the light bottom squark bi can be the lightest of all 
matter scalars. It has been shown that pairs of bottom squarks may be directly detectable 
at the Fermilab Tevatron pp collider if m^^ < 240 GeV, provided bi —>■ bZi is the dominant 
decay mode, m^^ is not too heavy, and sufficient integrated luminosity is obtained PB| , |53| . 
In Fig. 1^, we have denoted with stars in the lower left corner the points in parameter 
space where direct detection of bottom squark pairs should be possible, assuming 25 fb"^ 
of integrated luminosity at the Fermilab Tevatron. The region below the associated dashed 
contour shows where rrir < 240 GeV. 

01 



If /i > 0, for the allowed region shown in Fig. |TU|, bottom squarks are always too heavy 



to be detectable at the Tevatron. However, for some parameter space points, m^-^ < 210 
GeV, the putative reach of Tevatron upgrades [El]. For these cases, however, the associated 



value of m^^ from Fig. ^6 exceeds 140 GeV so that the top squark decay products are too 
soft to yield an observable signal at the Fermilab Tevatron [p^ ,^. 

We have also examined in a limited context the observability of signals in the RIMH 



model at the CERN LHC pp collider. In Ref. the reach of the CERN LHC for SUSY 
particles within the mSUGRA model has been examined. In that model, values of rrig ~ 
1500 — 2000 GeV are generally accessible assuming just 10 fb~^ of integrated luminosity. 
The strategy was to select events with 

• jet multiplicity, rijct > 2 (with -Erjet > 100 GeV), 

• transverse sphericity St > 0.2, 

• Er(ji), Et{j2) > and > E!r, 

where the parameter E^ is adjusted to optimize the signal. We classify the events by the 
multiplicity of isolated leptons, and in the case of dilepton events, we also distinguish between 
the opposite sign (OS) and same sign (SS) sample as these could have substantially different 
origins. For the leptons we require 

• prii) > 20 GeV (£ = e or /i) and MT{i,:^T) > 100 GeV for the li signal, and 

• Prih, ^2) > 20 GeV for = 2, 3, . . . lepton signals. We do not impose any Pt(^) > E^ 
requirement on the leptons. 



In Fig. |T7|, we show the maximum value of signal cross section (as) over square root of 



background cross section {^/aB) in fb^*^^, where the ratio is maximized over five choices of 
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= 100, 200, 300, 400 and 500 GeV. The numeral at each plot point corresponds to the 
value of which gives the maximum value, e.g. 1 corresponds to i?^ = 100 GeV, etc.. 
The SM background has been taken from Ref. [|26|. As before, the negative (positive) plot 



is made for a slice out of Fig. |^ (Fig. |ToD, where mi/2 = 0.25mi6 {mi/2 = 0.22mi6) along 



which S is large. The other parameters are shown on the figure. This plot can be used 
to obtain the SUSY reach of the LHC in each of the channels shown, for any integrated 
luminosity and any statistical significance of the signal. The dashed line labelled 10 fb~^ 



shows the 5a level for that integrated luminosity. In Fig. we have in addition required 
that as > 0.5 fb - if this is not the case for the E^ value that maximizes the statistical 
significance, we show a lower value of E!j^ in this plot. We have also checked that for the 
range of mig values for which the signal in any channel satisfies our "5cr" and "5 events" 
criteria, as/ as is automatically larger than 0.2. We see from frame a) that for the chosen 
slice of the /z < parameter plane, models mig ^ 2400 GeV should yield observable signals 
at the CERN LHC for just 10 fb~^ of integrated luminosity. This corresponds to a value of 
g ~ 1450 GeV, with nii^ ~ 450 GeV being the smallest of the squark masses. The best 



significance is achieved in the 0^ and 1£ channels. The corresponding results are shown in 
frame h) for the slice with /i > 0. In this case, mig ~ 2600 GeV can be probed for the 
same integrated luminosity, corresponding to uig ~ 1400 GeV, and the lighter top squark 
with mi^ ~ 625 GeV is the lightest squark. We emphasize that the cuts used above were 
designed to pick out signals from the cascade decays of heavy squarks and gluinos in the 
mSUGRA framework, and not for the IMH model with a relatively light bottom or top 
squark. Nevertheless, we felt that it is instructive to show these results because general 
purpose searches similar to these will almost certainly be the first to be carried out when 
LHC data become available. 

The reach can be somewhat extended if one concentrates on signals associated with direct 
top or bottom squark production. As an attempt in this direction, we require events with 
exactly two tagged 6-jets (with 6-tag efficiency of 50%), together with the following cuts: 

• > 100 GeV, 

• Pt(&- jet(l)) > 100 GeV, 

• pT{h-jet{2)) > 50 GeV, 

• + E EriQCD - jets) > 1500 GeV, 

where QCD jets must have px > 25 GeV. The SM background is assumed to come dom- 
inantly from ti production, for which we find aB{tt) = 1.26 fb. Our signal cross sections 
after these cuts are plotted in Fig. |TB|, for the same parameter plane slices as in Fig. |T^. In 
this case, for just 10 fb~^ of integrated luminosity, the reach in frame a) for /i < extends 
to mi6 ~ 2850 GeV, and in b) for yU > to mig ~ 2950 GeV. This corresponds to a reach 
in rrig ~ 1700 and 1600 GeV, respectively. We note that the rather hard cuts that we have 
made (to eliminate the top background) will likely cut out the direct ii or hi SUSY signal 
if the sbottom or stop mass is very small. In this case, dedicated searches for these may be 
possible. It should also be kept in mind that within the RIMH framework, t or 6 squarks 
may also be within the reach of future linear electron positron colliders. 
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VII. SUMMARY AND CONCLUSIONS 



An inverted scalar mass hierarchy has been suggested |TI|] as a possible way to amelio- 



rate the SUSY flavor and CP problems without destabilizing the weak scale. A particularly 
attractive approach pioneered in Ref. fl^ is to generate this hierarchy dynamically. These 
authors identified a simple set of SO{10) symmetric boundary conditions for the renormaliza- 
tion group evolution of the soft SUSY breaking parameters of an MSSM +RHN model, and 
showed that the large Yukawa couplings for third generation sparticles drove their masses to 
sub-TeV values, leaving first and second generation sparticle masses in the multi-TeV range. 
Our main goals in this paper were to examine the extent of the hierarchy that could be gen- 
erated within realistic RIMH supersymmetric models, and to explore the phenomenological 
consequences of these scenarios. 

In contrast to values of the crunch factor ( p.6| ) in the range S = 50 — 800 reported in 
Ref. [0, we found that a much more limited IMH could be generated, with S < 9. We 
found several reasons for the difference. First, in our bottom-up approach, starting from the 
observed values of matter fermion masses, to calculate third generation Yukawa couplings, 
only solutions with /{Mgut) ~ 0.5 could be obtained, to be compared with /{Mgut) ~ 1 — 2 
used in Ref. |T^. Using a top-down approach, we checked that such large values of the GUT 
scale Yukawa coupling lead to rrit ~ 195 — 210 GeV, which is well beyond its experimental 
value measured at Fermilab. We also verified that for /{Mgut) ~ 1-2, values of S up to 
35 could be obtained but, of course, with unacceptable values of m^. The second difference 
came from the inclusion of numerous perturbations to the simple analytical solution of Ref. 



14| which reduce the hierarchy that can be obtained. These include i.) SSB terms of order 



the weak scale in the ROE, ii.) splitting between Yukawa couplings and also A-parameters 
below the GUT scale, in.) two- loop contributions to renormalization group running, and iv.) 
inclusion of 5*0(10) D-term contributions to scalar masses: these contributions, which are 
essential [0] for obtaining solutions with REWSB, perturb the simple boundary conditions 
suggested in Ref. |1| . 



While these considerations make it clear why we are unable to find model solutions 
where third generation scalars are an order of magnitude lighter than their counterparts of 
the first two generations, we have found realistic models with REWSB that give first and 
second generation scalar masses of several TeV, at the same time maintaining sub-TeV third 
generation scalar masses. 

We attempted a systematic exposition of parameter space leading to RIMH models. 
The parameter space for which the largest hierarchy is obtained is characterized by large 
tan P ~ 50, Aq < 0, Md ~ 0.2mi6, and for any fixed value of mi/2, as large as possible values 
of mi6 as allowed by REWSB (see Figs. |^ and |TD|). The RIMH solution can be obtained for 
either sign of /z; solutions with negative /i tend to lead to a somewhat larger hierarchy. 

To analyze the phenomenological implications of RIMH models, we have selected two 
slices of the parameter space which yield large S, one for each sign of /z. The corresponding 
sparticle masses are shown as a function of mig in Fig. ^ for /i < 0, and in Fig. for > 0. 
The model predictions for the neutralino relic density are shown in Fig. |13[ We found that a 
significant fraction of the model parameter space is excluded because it leads to too young a 
universe. This is because sparticles are heavy so that the LSP annihilation cross sections are 
generally suppressed. The relic density prediction falls close to or within the experimentally 
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favoured region for models characterized by either a very hght bottom squark (yU < 0), or by 
a hghtest neutrahno in the transition region between being bino-hke and higgsino-hke (both 
signs of fj,). 

We have also evaluated the branching fraction for the decay 6 —> 57 as a function of mig 
for these same parameter space slices. Our result is shown in Fig. [1^. We see that models 
with negative /i yield too large a rate for this decay, and are essentially excluded by the 
CLEO experiment |]5D| at the 95%CL, except for the largest values of uiiq in the figure. We 
note, however, that the branching fraction is < 25% larger than this upper limit as long 
rriiQ > 4600 GeV. Since it is possible that other effects (squark mixing, SUSY phases) may 
modify the calculation of this branching ratio, it may be worthwhile to exercise some caution 
before definitively ruling out these scenarios. Models with positive /i can yield agreement 
with the experimental measurements for mie ~ 1800 — 2000 GeV. Unfortunately, this range 
is conclusively excluded by the relic density constraints, assuming that the LSP is absolutely 
stable (as it is within our framework). Combining the constraints from the 6 ^ 57 decay 
with those from the relic density, we see that values of niiQ ~ 4600 — 5000 GeV and /i < 
might lead to the most viable of these large S scenarios. 

Finally, in the previous section, we have examined some consequences of RIMH models 
for collider searches. For the Fermilab Tevatron, only very small regions of model parameter 
space could be accessed via direct searches for bottom squark pair production. At the CERN 
LHC collider, the limited IMH obtained means that a complete decoupling of heavier states 
does not take place, in contrast to what happens [1^ in GSIMH models. By searching for 

+ jets events with at least two tagged 6-jets, values of nig ~ 1600 — 1700 GeV could be 
probed with just 10 fb~^ of data. There are, however, potentially viable regions of parameter 
space (including the region favored by the combined constraint from the neutralino relic 
density and radiative b decay) where finding supersymmetric signals could pose a much 
bigger challenge. 
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Appendix: Two loop RGEs for the MSSM plus RHN model 

In this appendix, we give the relevant two loop RGEs for the MSSM plus right handed 
neutrino model adopted in this paper. We augment the field content of the MSSM by adding 
a singlet neutrino superfield for each generation i = 1, 2, 3. In keeping with the basic 
structure of 5*0(10) GUT models, we will assume that only the third generation neutrino 
Yukawa coupling is significant, and will neglect terms including neutrino Yukawa couplings 
for the first two generations. The form of the superpotential and SSB terms given in Eqs. 

and are then sufficient to determine the form of the renormalization group equations. 

The two loop RGEs for the MSSM+RHN model can be extracted from Ref. [Q. For 
the gauge couplings, we find 
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The expression for the gaugino masses can be written in terms of the same matrices of 
coefficients: 



dt 167r2 " (167r2)2 
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For the A-parameters, we find 



d . 
dt 



167r2 



(IGtt 



2\2 



with 



= 2 [E c^^^^M, + 6/,2A, + f^At + f^A^ 



P'a! = 2 [E + + + f^A^ 

where the Cj, and c'- are given in Ref. and c"' = (|, 3, 0). Also, 



/o(2) 



(2) 



225 



Ml (2743^7? + 225^72^ + + ml + 270/,^) gl - 1A, [fl + ft) fl 



-2glM^ (gl + Ibgl + ^gl + G/^) - ^^qIM^ [i^qI + 45(?2 - 20^73' + 90/^' 
+ (6^2 + 30^2 ^ gQ^2 _ 25 j2 _ 220/2 _ 15/2) 

+ \a,P, {2gl - 50fl - 25f^ - 5f^) - 2AJ'^ (sf? + Qfl + fl 
-^Mi (287^2 ^ 45^2 ^ 4Q^2 ^ ^g^2 ^ 3Q^2 ^ 54^2^) ^2 



-2(72M2 + i^gl + 8(71 + 6/2) - ^glM^ [gl + 9gl - Agl + 18/,^) 
+|a/,2 (4^? - 25/,2 - 50fl - 5fl) + ^Afr {Qgl - i5fl - 5/.^ - 30/; 



+-Ab/2 (2(72 + 30^72^ + 80^73' - 220/,^ - 25 fl - 15 fl) - 2AJl (fl + f 
Ml (135^7? + 9^72' - 2fl + 6fl) gl - 32glMsfl - 6AJl [fl + fl 



2' 

T I 5 



5 

+^A./2 [Qgl + 30(72^ - ^5fl - 15 fl - imfl) - 6AJl [fl + 2fl + fl 
-\glM^ {^gl + 25<72' + 10/,^) - \A,fl {2gl - mgl + 90/,^ + I5fl + A5fl 



--Ml (207^7^ + 45^72^ + 20// + 30/,^ + 30/,^ j gt - 32giM,fl 

+jAJl (Agl + SOgj - 15 fl - 90 fl - A5fl) - ^g^M, {3gl + 25gl + 10 fl 

+IaJI {6gl + 30gl - A5fl - 100 fl - 15 fl) + ^AJl {2gl - 5fl - 5fl 



-QAfl {fl + fl] 
For the soft SUSY breaking scalar masses, we find 
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where for the Higgs masses, we have 
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In the above equations the sums over color and flavor are included, and Tr simply means a 
sum over generations. The term 5 = if universality is imposed, but can be significant if 
scalar masses are non-universal. 

The third generation matter scalar /^-functions are given by 
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TABLES 



TABLE I. Weak scale sparticle masses and parameters (GeV) for two IMH model case studies. 
The first has /x < 0, and the second has > 0. 



parameter 


case 1 


case 2 


mi6 


2500.0 


2500.0 


mio 


3535.5 


3535.5 


Md 


500.0 


500.0 


Mn 


650.0 
1.0 X 10^ 


650.0 
1.0 X 10^ 


Ao 


-5000.0 


-5000.0 


tan/3 


50.0 


50.0 


rrig 


1547.4 


1558.9 




2789.4 


2784.2 




2589.9 


2585.5 




2367.7 


2366.1 




2566.7 


2567.1 




2366.3 


2364.8 




858.8 


742.0 




1250.6 


1500.4 




446.4 


1422.0 




1126.0 


1527.9 




1089.0 


661.2 




1299.0 


878.8 




1088.5 


862.4 




518.6 


528.0 


ruh 


518.5 
281.5 
129.1 


527.8 
281.2 
129.1 


rriA 


819.3 


1839.9 




826.3 


1844.2 




-651.5 


849.1 


R 


1.03 


1.74 


MMgut) 


0.535 


0.495 


MMgut) 


0.544 


0.305 


JAMgut) 


0.552 


0.531 


S 


7.0 


4.6 



27 



TABLE II. Weak scale sparticle masses and parameters (GeV) for an IMH model us- 
ing a top-down approach. In this case, the unified Yukawa coupling at the GUT scale is 
f(MarjT) = 1.282, and the top quark pole mass is found to be m,f = 204.8 GcV. We take /x < 0. 



parameter 


value 


mi6 


12766.3 


mio 


18054.3 


Md 


765.4 


Til 1 In 
'"'1/2 


732.3 


Mjv 


0.0 X iU 




-ZOOoZ. I 


tanp 


OZ. 1 


m-g 


ZUOZ.O 




1 9fi8/l 7 


TTlj 


lOROI Q 




1 ORRO K 
IZDDZ.O 


rrip 


lz779.z 








12662 3 


tl 


1 /1 1 Q n 


'"■t2 


1955.5 


1717 

bi 


1041.3 


mr 

'"^62 


1931.7 


'"'Tl 


3381.7 


Tn ~ 
1 1 i"f2 


3725 7 


•''■Ur 


^791 8 




735.7 






735.6 




376.4 


rrih 


153.6 


ruA 


779.0 


mH+ 


786.7 




-2367.7 


R 


1.0 


S 


33.9 
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Vi 





FIG. 2. A plot of crunch factor S versus a) rnxj^l'^XQ and 6) MnlmiQ, for all models with 
^0 < and /x < 0. 
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FIG. 4. A plot of crunch factor S versus a) rnxjilmx^ and 6) MD/miQ^ for all models with 
< and /x > 0. 
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FIG. 5. A plot of a) the physical value of nit versus /{Mgut), t>) crunch factor S versus 
/{Mgut), and c) S versus Mo/mie and d) S versus mi/2/mie, using a top-down solution to 
RGEs, and a unified GUT scale Yukawa coupling. In these plots, we adopt /x < 0. 
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FIG. 7. A plot of allowed regions for the RIMH model in the mie vs.mi/2 plane, for 
Mjv = 1 X 10^ GeV, Md = 0.2mi6, tan (3 = 50 and /U < 0. Points with triangles are excluded by 
a non-neutralino LSP, open circles by lack of REWSB and dots by tachyonic scalar masses. We 
show contours of 5 = 3 and 5, R = 1.05, rrig = 2 TeV, rriu^ = 3 TcV, fi = —1 TeV, and m^-^ and 
mjj^ of 1 TeV. The stars denote points where bottom squarks should be detectable at the Fermilab 
Tevatron, assuming an integrated luminosity of 25 fb~^. 
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FIG. 8. A plot of a) S, b) R, c) and d) m^^ versus mig for mi/2 = 0.25mi6, tan/? 
and /X < 0, for the various values (in GeV) exhibited in the figure. 
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FIG. 9. A plot of various sparticle and Higgs boson masses and —fx versus mie for 
mi/2 = 0.25mi6, tan/3 = 50, /x < and Mn = 1 x lO'^ GeV. 
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FIG. 10. Same as Fig. except ^ > 0. The degree of Yukawa unification ranges from 
R = 1.6 — 1.9 thoughout the plane. 
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FIG. 11. A plot of a) S, b) R, c) m^^ and d) m^^ versus ijiiq for mi/2 = 0.22mi6, tan/3 
and /X > 0, for the various Mjv values (in GeV) exhibited in the figure. 
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FIG. 12. A plot of various sparticle and Higgs boson masses and versus uiiq for 
mi/2 = 0.22mi6, tan/3 = 50, /x > and Mn = 1 x 10^ GeV. 
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FIG. 13. The neutralino relic density as a function of niiQ, for Mu = 0.2mi6, tan/3 = 50, 
Mn = 1 X 10^ GeV and a) /x < with mi/2 = 0.25mi6 and b) /j, > with mi/2 = 0.22mi6. 
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FIG. 14. The branching fraction for the decay b ^ for the same two sHces of the parameter 
plane as in Fig. 13 versus the parameter m-ig. The CLEO experiment requires that this branching 
fraction is between 2 x 10~^ and 4.5 x 10~^ at the 95% CL, while the recent measurement from 
the BELLE collaboration has a central value only slightly higher. 
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FIG. 15. The Wilson coefficients a) Cj, and b) Cs that enter the computation of the radiative 
decay of the b quark, versus the model parameter mie for the parameter space slice with negative 
fj,. The dashed lines labelled, SM, and W, respectively denote contributions to Cy or Cs from 
the SM tW loop, the tH~^ loop and the Wiijj loop. The sum of these is shown by the solid line 
labelled Ct^s{Mw)- The result of evolving this down to the scale mb is shown as the other solid 
line labelled Cj^si'^nb)- 
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FIG. 16. The same as Fig. except for the parameter space slice with /i > 0. 
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FIG. 17. A plot of signal cross section over square root of background cross section for 
Oi, li, OS and SS dilepton events after cuts from the RIMH model at the CERN LHC, versus 
the parameter mig. The n — a significance can be gained by multiplying by the square root of 
the integrated luminosity in fb. In a), we take mi/2 = 0.25mi6 and /x < 0. In &), we take 
mi/2 = 0.22mi6 and > 0. In both frames, we take Mo = 0.2mi6, tan/? = 50 and Mn = 1 x 10^ 
GeV. 
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FIG. 18. A plot of signal cross section for multijet + events which contain exactly two 
tagged 6-jcts, after cuts, from the RIMH model at the CERN LHC, versus the parameter mig- In 
a), we take mi/2 = 0.25mi6 and < 0. In b), we take mi/2 = 0.22mi6 and > 0. In both frames, 
we take Md = 0.2mi6, tan/3 = 50 and Mn = 1 x lO''' GeV. 
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